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Abstract 

Through an ^-expansion of the confined Calogero model with spin exchange interac- 
tions, we extract a generating function for the involutive conserved charges of the Frahm- 
Polychronakos spin chain. The resulting conservation laws possess the spin chain yangian 
symmetry, although they are not expressible in terms of these yangians. 
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1. Introduction 

Integrable spin chains with long-range interactions have remarkable properties, not the 
least being that they furnish a sort of discretization of particular conformal field theories 
with Lie group symmetry [1, 2, 3]. The archetypal model is the Haldane-Shastry model 
[4] in which N su(n) spins placed equidistantly on a circle are coupled by a spin-exchange 
interaction proportional to the inverse square of their chord distance: 

i N 

Here, Zj = exp (i^jf-), Zij = Zi — zj and Pij is the operator which exchanges the i th and j th 
spins. The primed sum indicates that the summation variables are restricted to differing 
values. 

The Haldane-Shastry model possesses a yangian symmetry algebra which can be taken 
as a manifestation of its integrability [1] . The conserved charges directly associated to this 
symmetry are not scalar (they transform in the fundamental representation of su(n)) and 
do not commute among themselves (they generate the yangian algebra, which is non- 
Abelian). However, from these charges, one can build a set of N scalar commuting opera- 
tors which turns out to be directly related to those obtained in [5] . However, this set does 
not explicitly contain the Hamiltonian, contrary to the natural expectation. Moreover, 
two additional conservation laws were known from brute force calculations [1, 6] but did 
not appear in this sequence. One expects that, together with H^ ns \ these represent the 
first few of a new sequence of a conserved charges. It is natural to try to fit this other 
sequence in a general scheme based on the fundamental object at the root of integrability 
: the monodromy matrix. For the Haldane-Shastry model, this has been accomplished by 
Haldane and Talstra [7]. They showed that the 'new' conservation laws can be obtained 
by taking a rather subtle limit of the more general dynamical spin model. 

For the well-known XXX model, which has short-range interactions, there are also two 
sets of conservation laws: there is a yangian symmetry [8] , out of which scalar conservation 
laws can be constructed and, in addition, there is a sequence of conservation laws that 
includes the Hamiltonian [9] . These two types of conservation laws are easily distinguished 
in models with short-range interactions: the first set is non-local (i.e., the conserved charges 
involve interactions of all the spins and they become truly non-local in the continuum limit) , 
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while the set containing the Hamiltonian is local (i.e., the n-th member of this sequence 
has a leading term describing the interaction of n adjacent sites). 

For spin chains with long-range interactions, the distinction between locality and non- 
locality is rather artificial, both sets of charges being manifestly non-local. The difference 
between these two sets lies in the fact that the Hamiltonian set found by Haldane and 
Talstra commutes with the symmetry algebra while the yangian set does not. Since both 
sets commute and can therefore be simultaneously diagonalized, this means that the eigen- 
values of the Hamiltonian set are degenerate and characterize a given multiplet while those 
of the yangian set can be used to label the differing states inside the multiplet. 

Let us point out, en passant, another major difference between integrable long- and 
short-range interacting chains, apart from the relativity of the locality concept. For short- 
range interacting chains, there exists a boost operator that allows for a recursive construc- 
tion of the local conservation laws. Its origin can actually be traced back to the transfer 
matrix formalism and the locality of the interaction [10]. No such operator is known for 
long-range interacting chains. 

The argument of [7] relies on a limiting formulation of the Haldane-Shastry spin chain. 
The model can be viewed as a special reduction of a general Sutherland model (a dynamical 
model with sin -2 r interaction) with spin degrees of freedom. 

The introduction of the spin degrees of freedom in a Calogero-Moser-Sutherland model 
is rather direct [11] (see also [12]) . If in the classical version of the model, the potential 
takes the form gf(ri, rj) (up to a possible harmonic part), where g is a coupling constant, 
the quantum version reads 9(9 + l)/( r i; r j)- The integrability turns out to be preserved 
if the term g(g + 1) is replaced by g(g + Kij) where interchanges the positions i and j. 
The spin degrees of freedom can be introduced directly by imposing the K t j to be a spin- 
exchange instead of a position-exchange operator. Another approach, albeit less direct, 
amounts to retain the position meaning of but consider states that are symmetric 
under the interchange of both the position and the spin variables. The resulting effect is 
identical. 

The transition from a dynamical model with spin degrees of freedom to the spin chain 
has been phrased in general terms by Polychronakos in [13]. The idea is simply that from 
a dynamical model with spin degrees of freedom, we can somehow freeze the latter to 
generate a spin chain. However, this freezing entails a compatibility condition that follows 
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from the original equations of motion: the position variables must correspond to the zeroes 
of the potential. For the sin -2 r interaction potential, this fixes the positions of the chain 
sites to the roots of unity. Note, on the other hand, that if the potential contains an 
harmonic piece, this part does not contribute to the spin interaction potential but it enters 
in the definition of the minima (in fact, whenever it is present, it ensures the existence of 
these minima). 

In this letter, we study the Hamiltonian conservation laws of the Frahm-Polychronakos 
spin chain [13, 14]. It originates from a Calogero model with inverse square interaction and 
an harmonic confining potential, augmented with spin degrees of freedom. The potential 
minima fix the sites of the chain to correspond to the zeroes of the Hermite polynomial 
Hn(x), to be denoted Xi. The Hamiltonian takes the form 

i N i 

^ i j — \ ^IJ^Jl 

and we will consider the general case of su(n) spins, each of the iV spins belonging to the 
fundamental representation. This model has already been shown to be integrable and to 
possess a yangian set of commuting operators [13]. In the following, we will show that 
Haldane-Talstra's argument, formulated here in a somewhat different way, can also be suc- 
cessfully applied to this model, effectively generating the set of Hamiltonian conservation 
laws. 

2. Integrability and Conservation laws 

2.1. The yangian algebra Y[su{n)] 

Let us first briefly review the Yangian algebra Y[su(n)] (see for instance [8, 15]) focus- 
ing on its relation to the construction of commuting invariants. For all known integrable 
spin chains (except, in fact, for a single and somewhat pathological example [16]), the in- 
tegrability property can be traced back to the existence of a monodromy matrix, annxn 
matrix of operator entries which depends on a spectral parameter u and which satisfies 
the RTT relation: 

R(u - v)tW(u)tW(v) = tW(v)tW(u)R(u - v) . (2.1) 
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Here, the superscripts refer to two auxiliary subspaces in which the matrices act non- 
trivially, e.g. 

TW(tt) = T(«)®l nxn (2.2) 

and R, called the R-matrix, is an n 2 x n 2 c-number matrix which must satisfy the quantum 
Yang-Baxter equation: 

R( 12 \u)R( 13 \u + v)TL( 2S \v) = R( 23 )(t>)R (13) (« + v)R( 12 \u) . (2.3) 

The RTT relation ensures that the transfer matrix t(u), which is defined as the trace of 
the monodromy matrix t(u) = Yla=i T aa (u), satisfies 

[t(u) J t(u)] = (2.4) 

so that its expansion in power series in u~ l generates commuting conserved quantities (see 
below) . The Yang-Baxter relation is simply a compatibility relation for the RTT relation. 

A simple solution to the Yang-Baxter equation is given by Yang's rational solution 

R«J)(u) =u + AP (ij) , (2.5) 

where A is an unspecified deformation parameter and pW) exchanges the auxiliary sub- 
spaces i and j 

p(12) A (l) B (2) = B (l) A (2)p(12) _ (2 _ 6) 

With this choice of R-matrix and with the monodromy matrix expanded in a Laurent 
series as (denoting the ab matrix entry of T(u) as T ab ) 

oo 

T ab {u) = 5 ab + A u~ {rn+1) T^ , (2.7) 

m=0 

the RTT relation reduces to the following commutation relation 

i-\ 

\ rpab rpcd~\ sadrpcb rcbrpad , \ \ ^ J rpcb rpad rpcb rpad 1 (n a\ 

[ 1 £ ) L m J — 1 e+m~° 1 e+m + A 2_^\ k+m 1 £-k-l ~ 1 e-k-l 1 k+mj ■ 

fc=0 

From this structure, we can define two sets of commuting operators. One of these is 
obtained by the spectral expansion of the transfer matrix 

n 

[J*,J m ]=0 , I m = Y,T%. (2.9) 
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The other set is related to the quantum determinant of the monodromy matrix [17] 
Det q [T(w)] = J2 <o)T la{l \u-{n-l)\)T 2a V\u-{n-2)\\gr)...T n ^ n) {u). (2.10) 



crESr, 



Here, a(i) is the image of % under the permutation a, e(a) is the permutation's parity and 
the sum is taken over all permutations of (1 ... n). The quantum determinant is analogous 
to the Casimir operator of a Lie algebra in that it commutes with all generators: 



Det q [T(u)] , T(v) 



0. 



(2.11) 



This property allows one to define a second set of commuting operators from the coefficients 
of the series expansion of Det q [T(w)] in terms of the spectral parameter 



[j m ,J £ ]=0 , Det q [T(u)]=l+J2 u ~ {m+1)j m- 



(2.12) 



m=0 



A given Hamiltonian H will therefore be shown to be integrable if one can prove its 
symmetry under a non-trivial monodromy matrix 



[H, T(u)] =0 



(2.13) 



which guarantees the conserved character of the involutive sets I m and J m . However, 
monodromy matrices are formidable objects which usually do not allow such commutators 
to be directly carried out. It is therefore very useful to codify the monodromy matrix 
in a minimal form. Such a minimal coding can sometimes be realized in terms of the 
yangian algebra y[su(n)]. In fact, defining the lower-order yangian generators (Qo b ,Qi b ) 
(a, b = 1 . . . n) by i 

A 



Qo 



rpab 

~ J 



Qi 



ab 



-Tf+--(T T ) . 



\ ab 



(2.14) 



the first few of the commutation relations (2.8) read 

[Qo 6 , Qo d ] = S bc Qo d - S da Qo 
[Qo 6 , Qf } = S bc Qf - 5 da Q? 
Qo , [Qi , Qi ] — Qi , [Qo , Qx ] 

A " ([ Qo b , [(QoQo) cd , (QoQo) e/ ]l - (QoQo) a& , [Qo d , (QoQo) e/ ] 



;2.15) 



Here and hereafter, we use the obvious matrix notation (ToTo) ab = ^ ra _ T TQ C T§ b 
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These three relations define, or more precisely, completely characterize the yangian algebra 
Y[su(n)]. The third relation is a sort of compatibility requirement on the different ways 
to reach Q 2 from multiple commutations involving lower-order charges. 

One can reconstruct the whole monodromy matrix strictly from its lower-order yan- 
gians whenever the former possesses a trivial quantum determinant 

Det q [T(u)] = c-number . (2.16) 

To justify the last statement, consider the following special cases of the algebra (2.8) : 

= [ , TT ] + A(T-T od - T?Tg) (a ? d) 

7T +1 - T- +1 = [T- , T? c ] + A(T-T aa - T CC TD (no sum) . ^ 

The first of these relations allows us to compute any Tfi~{a 7^ b) in terms of the lower- 
order generators but the second relation is not sufficient to compute the T^ a by recurrence. 
However, if (2.16) holds, its spectral expansion gives a set of conditions on Y^a=i ^C" 
which, when supplemented by (2.17), allow one to compute all the from the lower- 
order yangians. This allows for a tremendous simplification because the symmetry of H 
under T(u) then follows from its symmetry under the induced y[su(n)] representation 

if{Det q [T(u)] = c-number j then j [H , Q$ 1} ] =0 [H , T(u) ] = o} . (2.18) 

When dealing with an irreducible representation of T(w), the quantum determinant must 
necessarily be proportional to the identity and the monodromy matrix can then be rep- 
resented by its lower-order yangians. However, when the considered representation is 
reducible, one must exercise care because the quantum determinant may then be a non- 
trivial operator. In the case of the Frahm-Polychronakos model, we will see that the 
symmetry algebra is reducible but nevertheless possesses a trivial quantum determinant 
so that in this special case (and for the Haldane-Shastry model), the monodromy matrix 
will be solely expressed in terms of its reducible lower-order yangians. 

2.2. The yangian representation in terms of Dunkl operators 

Having discussed the theory of su{n) yangians, we now focus on the construction of 
specific realizations useful for long-range interaction models. First of all, we work in a 
Hilbert space of iV particles endowed with su(n) spin, in which the position (momentum) 
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operator of particle % will be denoted by Ri (Pi); its spin operators are chosen to be the 
n 2 fundamental generators Ef 3 (a,b = 1 ... n) satisfying 

[ Ef , Ef ] = 5ij (5 hc Ef - 5 ad Ef) . (2.19) 

We now define an hermitian exchange operator K^, which permutes the positions of 
particles % and j 

k tj | r^... rfl. rf... r^ > = | rfl rf... rf... > . (2.20) 

We stress that in our notation, operator subscripts refer to particles whereas ket subscripts 
(superscripts) refer to positions (particles) so that 

Ri | rf... rf... r? > = r p | rf... rf... > . (2.21) 

The permutation operator Kij satisfies 

kijfiiu, Pi) = f(Rj, Pj)kij k l3 k Jk = k lk k l0 

^ £ _^ ^ ^ ^2 22) 

k ij f(R e ,p e ) = f(R e ,p e )k ij k ij k ki = k ke k i j 

Here, the caret is used to stress that k^ is an abstract Hilbert-space operator and therefore 
acts trivially on any c-number. This contrasts with the operator generally used, which 
exchanges the position eigenvalues according to Kijri = rjKij and which is simply the 
position-space representation of this abstract operator: 

<r[ 1 Lr^ ) \k ij \^>=K i j<r[ 1 Lr ( ^ ) \^> . (2.23) 
In a similar way, one can also define a spin exchange operator 

p Tpab pa6 p 

, , • (2-24) 

Pi 77'ao 7-rao r> 

In the fundamental basis, this operator takes the simple form 

n 

Pij = E t E f • (2-25) 

a, 6=1 

Now, in order to eventually establish a link between spatial and spin models, one 
introduces a projection IT [18], which consists in projecting onto states that are symmetric 
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with respect to the joint interchange of position and spin variables, that is, states satisfying 
KijPij = 1. In practice, this projection boils down to the following operation: in a given 
expression, we move all operators to the right and replace them by Pij operators acting 
in reverse order; e.g. 

n {kijkjk} = p Jk p tJ . (2.26) 



This projection possesses the following crucial properties 

U{AB} = U{A} ■ U{B} if [ B , K^Py ] = 

U{AB} = U{B} ■ U{A} if [A,U{B}]=0. 



(2.27) 



Using this projection technique, one can construct a spin representation of the algebra 
(2.8) [19]. This representation is based on given position-space Dunkl operators Di (and 
from now on, we will use the overhead bar to indicate that an operator acts non trivially 
only in position space) obeying 



[Di,Dj] =X(D i -D j )k ij . 
By induction, one can in fact prove the more general commutation relation 



(2.28) 



i— i 

[Df ,Df]= A^( J Df +m J Dj- fc - 1 - D t - k - 1 D'< +m ^k ij , (2.29) 



k=0 

from which one can immediately define the following involutive set 

N 



[l e ,lm]=0 , I m = Y, A m • (2-30) 



i=l 



These quantities are purely spatial; in order to define spin invariants, we can use the 
properties (2.27) and (2.29) to show that the currents 



N 



T% = Y,Et a Tl{D?} (2.31) 



i=l 



satisfy the monodromy matrix algebra (2.8). The involutive I m set associated with this 
algebra is then simply given by 



N 



[l m ,I t ]=0 , I m = U{I m } = J2^{DT}- (2-32) 



Quite remarkably, the monodromy matrix (2.31) can also be expressed in the form T(u) = 
U{T'(u)}, with T'(u) another representation of the algebra (2.8), given by [18] 



1 N N 

T ' W = AMn{( n -^ + AE ' A^^Hiu-D',) (2.33) 



i=l 



where the spin operators have been grouped in matrix form according to (Ej) ab = E\ a 
and the product over % includes both the usual linear matrix product as well as the tensor 
one. Here, the D[ are modified Dunkl operators 



N 



Dl = A - A k H . ( 2 - 34 ) 

j<i 

which satisfy the degenerate affine Hecke algebra with respect to position-space permuta- 
tions 

Ku±iD' i - D'^Ka^x = ±A 

(2.35) 

[D>,D>}=0. 
Using this Hecke algebra, one can show that 

[ A(u) , K ii+1 } = [ A(u) , Kii+iPu+i } = , (2.36) 

which can in turn be used to prove that[T'(w) , K ii+ iPii + i ] = 0. But since any permu- 
tation can be expressed as a product of transpositions, we actually have 

[ A(u) , Kij ] = [ A(u) , KijPij ] = [T'(u) , KijPij } = . (2.37) 

Now the quantum determinant of T'(u) having already been calculated as 

Det q [T'H]=^±^, (2.38) 

one can use the property (2.37) to factorize the projection and calculate the quantum 
determinant of T(u) as [18] 

Det q [T(«)] = Det q [n{1»}] = n{Det q [T'(u)] } = n| A ^ ( + A) | . (2.39) 



In principle, one can extract the {Jk} set from this formula but the result is highly cum- 
bersome. It is much simpler to focus instead on A(u). Indeed, the relations (2.35) and 
(2.37) also imply 

[A(u), A(v)} = [A(u), A(v)} =0, (2.40) 

where A(u) = n{A(-u)}. One can therefore define a simple set of commuting operators by 
using Jjj In [A (it)] as their generating function [7] 

N 

[H m ,H £ }=0 , H m = J2{Dd m 

1=1 (2.41) 

[H m ,H t ]=0 , H m = Tl{H rn }=Y,n{{D' l ) m y 

i=i 

By virtue of (2.39), this new set is obviously equivalent to the {J m } set and from now on, 
we will focus on the sets {I m } and {H m }. 



3. The Dynamical Calogero Model 

In order to obtain the Hamiltonian spin-chain conservation laws, one must first con- 
sider an iV-body dynamical Calogero model in which the particles are chosen to have 
unit mass and are allowed to move along the line, under the influence of a position-space 
exchange interaction and subject to an harmonic confinement 

1 N i N ( £ ir \ i N 

h (CC> ^ \ E p ' - \ E ( „ p + & E *? • p-d 

z i=i i,j=i z i= i 

The integrability of this model has been demonstrated e.g., in [20] by means of the operators 

N 1 

f = Pj + V^ig JT k i k ± ^^ r j > (3-2) 

fc=i ^' fc 



X>7 



which satisfy (2?- ) = £>J, in addition to the commutation properties 
[t>f,Vf] =0 

[Pf , Pf ] = T 2w(Jy (ft + <? £ £ ifc ) ± (1 - 5 l3 )2ugK 

^ fc=i ' 



(3.3) 
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From these, one can define [21] the Dunkl operator Pj = T>^~V i + Tllo 
N 1 N f 1 1 

v * = ~ g2 E inr + E ( p * + p ti - ug r 

^! ^Ik-tiki ^ I J 



J,K=L ^ J = 



3=£i 



(where here and hereafter, we use the notation if^...^ = n?=i ^jij+i)- The deformation 
parameter of this Dunkl operator is A = —2ujg; it therefore satisfies the commutation 
relation: 

[Pi , Vj ] = -2ug(T>i - Vj)kij . (3.5) 
The representation of the yangian algebra induced by this Dunkl operator is given by 

N 



i=i 

N N , ., r j s 

<W yj' ( E,E 3 E t )^ + y:'W + (36) 

i,j,k=l J J i,j=l J 

i=i I j=i y J 

Its associated monodromy matrix then allows us to generate two non-trivial involutive 
sets of operators, denoted X m and H m (calligraphic symbols being used for the charges 
pertaining to the dynamical model). Their first member is given explicitly by 

JV N 

ii = E^ = 2 ^ (cc) - ug E k *j ( 3 - 7 ) 

i=l i,j = l 

N , N x 

^ = E ( ^ + 2 ^ E ^ ) = 27{{cc) ■ ( 3 - 8 ) 
i=i ^ j=i ^ 

Because H^ cc ^ is included in {H m } and these essentially arise from an expansion of the 
quantum determinant, the symmetry of H^ cc ^ under the monodromy matrix induced by 
this Dunkl operator is manifest. This means that both the T m and H m sets constitute 
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involutive invariants for this dynamical model. Moreover, the basic relations (3.3) can be 
shown [20] to imply 



thereby furnishing a set of creation operators from which the spectrum can be readily 
obtained. 

4. The Frahm-Polychronakos spin chain 

The Frahm-Polychronakos spin-chain model is defined by (1.2). In this expression, 
the XiS are the zeroes of the Hermite polynomials. In this section, we will consider a 
generalized version of the FP Hamiltonian H(r): 



in which r is a set of unconstrained position eigenvalues, in order to see explicitly how the 
yangian symmetry picks up the particular FP model, i.e. how it enforces = a^. In the 
following, the (potential) conserved charges pertaining to this general version of the spin 
chain that are inherited from the dynamical model will be denoted by Im'(r) and Hm\r). 
The subindex refers to an ^-expansion to be explained shortly. 

In order to generate a candidate symmetry algebra for this generalized model, we 
consider the position-space representation of the dynamical Calogero model, in which, 
from now on, we set u = g = 1. The spin part of 7i( cc ) is simply isolated as the linear 
h-piece of the Hamiltonian. More generally, the spin part is obtained by differentiating the 
Hamiltonian with respect to % and, since the kinetic term is quadratic in h, setting h = 
at the end. This ignores the fact that the zeroes should be fixed at particular positions, but 
nevertheless suggests to consider the ^-expansion of the Dunkl operators and the related 
conserved operators. 

Let us then expand the dynamical Dunkl operator (3.4) according to T>i — T>^Ti k 
and for the time being, concentrate on the zeroth order term: 




(3.9) 




(4.1) 




(4.2) 
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Since the Dunkl algebra (3.5) is satisfied for all values of h, V^^r) is also a genuine Dunkl 
operator, with deformation parameter A = —2. Now for generic values of the r^'s, the 
induced y[su(n)] representation is irreducible and its quantum determinant is therefore 
a trivial c-number. As a corollary, the {Hm\r)} do not provide non-trivial conserved 
charges, i.e. these quantities are independent of any exchange operators. On the other 
hand, the set {Tn\r)} does provide a non-trivial involutive ensemble, its first member 
being given by 

N N N 1 JV 

Z<°>M = E n H 0) (r)} = - £' P« + E' k + • 
To obtain this result, we used the identity 

^/ i i ^ r i i ii 

V = -.y { + + 1=0. 4.4 

jM=i r i kTkt 3 j^ =1 l r Jk r ke r ke r ej r ej r jk J 

Now, the higher order Im (r) do not contain any term having the form of H(r). In other 
words, this set of commuting operators has no relation at this point with the generalized 
model defined by the Hamiltonian H(r). In order for the set {I^(r)} to represent involu- 
tive invariants for (4.1), one must enforce the invariance of H(r) under the corresponding 
y[su(n)] algebra, whose lower-order generators are given by 



eg»(r) = -y> 



N 

T 

/ j — i 
i=l 

N N ( N -\ (4.5) 

i,j,k=i Tljrjk ~ 

In other words, we require that [H(r) , i)( r )] = 0. A direct calculation [21] shows 
that this holds if and only if the variables r; obey 




E i = h ■ (4.6) 



r 



One can show that this condition is satisfied by the zeroes (written Xi) of the Hermite 
polynomial Hn(x) (cf. Appendix B). In fact, by judiciously substracting known summation 
identities [22] for these numbers, one can generate a whole sequence of 'higher order' 
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identities, the simplest of them being listed in Appendix B; these will play a crucial role 
in subsequent calculations. 

In retrospect, by freezing the positions of the particles on the zeroes of H^{x), we 
send H(y) on _ff"( FP ) and thus obtain an integrable Y [sw(n)]-symmetric spin chain with a 
non-trivial involutive set of invariants given by 

N ( } 
7 m (x) ee lim 2g)(r) = £n|(£< 0) (x)) j . (4.7) 

These are the conserved quantities first found by Polychronakos [13] (but without the 
yangian interpretation). 

Moreover, defining ee Ym=i ^iipt)™ > expanding (3.9) to 0(%) and setting u = 
g = 1, we find 

[tf<cc)<o) (r) t c ±(i) (r)] + [tf(cc)(i) (r) j c ±(o )(r)] = ±mC ±(°) (r) . (4.8) 

Because ?i^ cc ^ ^(r) is scalar and X^-invariant, the first commutator on the left hand side 
reduces to the action of the derivative on H^ cc ^°\r), which is given by 



This vanishes asr^x (cf. the identity (B.7)) and (4.8) takes the form 

i N i 

[^e'tV^-^w] =±mC - 0)(x) - (4 - io) 

Taking now the projection and using the KjjP^-invariance of the two commuted operators, 
we obtain a whole set of creation operators for the FP model: 

( N N , 

[# (FP) ,C±]=±mC± , C±eeH J> fc £ — K M ±x k ) . (4.11) 

U=i e=i ^ kl ' ) 

These generalize the lower-order creation operators found in [21,14]. We therefore possess a 
set of non-trivial creation operators C^j and conservation laws I m . However, as previously 
pointed out, the {Hm} set associated with the symmetry algebra is trivial. Since {1$} 
does not contain the defining Hamiltonian, a whole set of commuting conservation laws is 
still missing. 



N ~\ N 
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5. The Hamiltonian conservation laws of the FP model 

Our proof for the commutativity of the conservation laws will strongly rely on the 
structure of the FP Hilbert space. For the su(2) Haldane-Shastry model, the yangian 
symmetry algebra has been shown to be a direct sum of irreducible y[su(2)] "motif 
representations, each possible motif appearing with unit multiplicity [18] . This result has 
been obtained by calculating the dimensions of the F[su(2)] motif representations as a 
tensor product of su(2) spin representations and then showing that these motifs exhaust 
the Hilbert space. For the su(n) case (n > 2), the motifs are not expressible as a free tensor 
product [21] and to our knowledge, it hasn't been proved that the Y"[s , u(n)] motifs exhaust 
the Hilbert space. However, strong numerical evidence [21] suggests that the symmetry 
algebras for both the su(n) HS and FP models are also a direct sum of non — degenerate 
motifs. In the following, we will consider this statement to be true. 

The non-degenerate character of the motifs implies that any two operators A and B 
commuting with the monodromy matrix T(u) must also commute amongst themselves (see 
e.g., [7]). Indeed, the Hilbert space of our reducible Y~[su(n)] invariant theory contains a 
certain number of yangian highest-weight states, each of which is associated with a given 
motif. These highest- weight states are eigenvectors of the diagonal elements T aa (u) (a = 
1 . . .n), with eigenvalues that completely specify the given motif. Since the motifs have 
unit multiplicity, the highest- weight states T aa (w)-eigenvalues form non-degenerate sets. 
Now consider the two states AB | A > and BA | A > , where | A > is a yangian highest- 
weight state. Since A and B commute with T(it), both these states will be eigenvectors 
of T aa (u) with the same eigenvalue. But since these eigenvalues are non-degenerate, the 
two states must in fact be proportional to one another, which implies [A, B~\ = on any 
highest- weight state. In fact, since all of the states can be generated by acting on the 
highest-weight states with lowering operators of the form Y[iT aibi (K) (with < hi and 
the Xi chosen to satisfy a set of Bethe ansatz equations), one sees that A and B will in 
fact commute in the entire Hilbert space. 

We will now prove that the first order terms in the /i-expansion of the dynamical 
{H m } set satisfy the FP Y[su(n)] symmetry and are therefore in involution. Starting from 
the dynamical symmetry [ <2(o,i)(r) , H m (r) ] = and expanding to 0{Ti), we have 

> = -[Qg 1} (r) , n^(r)] . (5.1) 
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On the right hand side, the commutation with Qq 1 "* is trivially zero, while that with 
can be greatly simplified by appealing to the scalar nature of Hm\r) and its invariance 
under K i3 - and K^P^f 

[Qf\r),«W(r)]=f'(E,E,) T ; l{A(w<°'W)+^(H!,° , W)}. (5.2) 

i,j — 1 

To further simplify the right hand side, we will now explicitly calculate (r). To this end, 
let us return to the abstract Hilbert-space formalism and consider the following integral 

F m (r) = Jd yi ...dy N < yi 1 }.. y^ | H^(Ri . . . R N ) | sym(r) > , (5.3) 

where | sym(r) > = J2 creSN \ r ^i) ■ ■ ■ r i^/v) > ano - the notation Tim (R\ . . . Rn) is used to 
stress that these charges are K i3 -independent. Applying Tim(Ri ■ ■ ■ Rn) to the left yields 



l)--- r a(N)> 



F m(r)= £ J dy 1 ...d yN Ht\vi...yN)<y { fi..y { N ) \r% ) 

= J dyi ■ ■ -dyNKmiyi ■ ■ -yN)8(yi - r a{1) ) . . .8(y N - r a(N) ) (5.4) 
= E ^in ) (r a(1) ...r a(N) ) = N\n^(r 1 ...r N ), 

ct&Sn 

where we have used the Ki 3 -invariance of Tim (r± . . . r^) in the last step. We therefore 
have the following result 

n^(r) = ^F m (r). (5.5) 
On the other hand, going back to (5.3), one can express F m (r) in the form 

/N 
d Vl ...dy N <y{ 1 l.y N V) \ ^{v'™}™ |sym(r)> , (5.6) 
i=i 

where the modified Dunkl operator has the following explicit expression 

JV N N N 

p/(0) = v (0) + 2 J- k l3 = E -f^~K ijk + E sgn(* - j)K tJ + £ w + R\ . (5.7) 

3=1 J,k=l J 3 = 1 3 = 1 l J 

Applying now Tim (Ri ■ ■ ■ Rn) to the right and using 

| sym(r) > = | sym(r) > (5.8) 

i) r W 

r(l)--- r a(N) 



f(Ri) |sym(r)> = £ f(r a(i) ) \ r% . . . > , (5.9) 
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one obtains 

N ( N 



w=EEI E 



ctESjv i=l v j,k=l 



ra(j)a(k)r a (k)a(i) 



N N ^ m 

+ + — • 

7 = 1 1 = 1 O-W^O) J 



(5.10) 



j=i j 



Considering now the (iV — 1)! permutations for which a(i) = £, this can be rewritten as 

N N ( N N N 1 ~] m 

Finally, using a binomial expansion to factorize the sgn(i— j) term (i.e., the r^-independent 
piece) and using (5.5), we find 

N m ( N N }P 

«8>m - E E cj E ^ + E ^ + . (5-12) 

fcip=o I j,fc=i Jfc fe£ j=i ^ J 

where 

Cp ™ = iv \v) ^i 2t ~ {N + l) T~ V • (5 - 13) 

To complete the calculation of the commutator (5.1), we need to evaluate the action 
of the derivative on Hm\r) (cf. (5.2)) and freeze the particle positions: 



F) ( -\ N m ( N _1 N 1 >[P ~ 1 

Bm£ {«£)« =jsEE^ E ^ + E^+*? 

2 k ^ £=i P =o I i,/c=i JK Kt i=i j 



j,fc=i j " ■■ i=i ^ 



^ I J=l r i* r " + ^ r ;< +r ) 



(5.14) 



In such calculations, the r — > x limit is not a simple substitution and must be taken with 
care. Indeed, the summation formulae for X{ are valid for numbers and are therefore not 
preserved by the action of the derivatives. This means that one may take the substitution 
r — > x and use the simplifying identities only if the targeted expression is no longer acted 
upon by any derivatives. In light of this remark, we see that we cannot simplify the second 
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bracketed factor in (5.14) without first carrying out the differentiation. However, we can 
use the formulae (B.7), (B.8) and (B.ll) to reduce the first bracketed factor right away: 

N , N x JV N 

E— ( - E— + — J + E 4- +*? = - £— +*? 

„2 



= -xj + {N - 1) + xj = N - 1 . 



We are thus left with 



N m r, ( N N ~\ 



(5.16) 



Commuting the sum over £ past the derivative and using the identity (4.4) we finally obtain 
d 



lim 



r^x dr 



n m p. ( N 1 AT ^ 

"J? w = E c - ^ - ir 1 ^ E ^ + E ^ + 



J2C pmP (N -1)^2 \x i -2j2-^ =0 

p=0 ^ 7 = 1 *J J 



(5-17) 



where we have used relation (B.9) in the last step. Substituting this result back in (5.2), 
we see that the quantities H 2m = lim r ^ x Urn (r) satisfy the FP Y[su(n)] symmetry. In 
other words, the i?2m all commute with the monodromy matrix. As already pointed out, 
this implies that they are necessarily in involution. 

A compact but implicit expression for the i?2m following from the ^-expansion of the 
dynamical operators, is given by 

N m—1 , ^ 

H *™ = Mm £ n Al{v) Bi {r)Ar p -\r) , (5.18) 
r x i=i p =o < > 



where 



Mr) = V —K ijk + V sgn(i - j)Kij + V \- + r, 

j,k=l Jk kt j=l j=l ^ 



(5.19) 



r- 
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The first two members of this set can be calculated as 

H 2 = V' PJ (5.20) 

iJ=1 x *3 X Jt 

Ha = ^2 \ (Pijkl 

i,j=i ^ i j l j i j j 

Because H 2 = 2H^ FP \ these involutive quantities are necessarily invariants of the FP 
model. Note that in the calculation of H 4: we have used the summation identities (B.8), 
(B.15) and (B.16) in order to simplify the final result. We have also been able to crosscheck 
the commutativity of H 2 and H4 by a direct computation, the details of which are given 
in Appendix A. It should also be pointed out that since it is impossible in quantum 
mechanics to freeze the particle positions onto the lattice sites and enforce at the same 
time the vanishing of their momenta, the absence of derivatives in H 2 and H4 should 
not be regarded as a consequence of the freezing procedure but as a rather impressive 
mathematical cancelation, whose raison d'etre has yet to be determined. Also, notice 
that our limiting procedure does not generate any odd-type conservation laws because the 
dynamical Calogero model simply does not possess such symmetries. On the other hand, 
it is easy to see that the FP model does possess such symmetries by verifying explicitly 
that it commutes with the following operators 

* = < 6 - 22 > 

i=i 1 

^±{^)^-\£(^)l- (-3) 

Here, we have not performed any simplifications on H3, in order to show that [Hi , H% ] = 
(recall that the summation identities can only be used once all derivatives have been 
commuted to the right). It therefore seems possible to generate odd- type H m , albeit by 
brute force. These odd conservation laws seem to commute amongst themselves as well as 
with the even H m although they manifestly do not possess the yangian symmetry. This 
means that Haldane and Talstra's argument cannot be used to isolate their generating 
function; we have not yet found the generating function for such a set. 
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6. Conclusion 

Using an ^-expansion of the dynamical Calogero model, we have succeeded in con- 
structing an even set {H 2 , H 4 . . .} of involutive charges for the Prahm-Polychronakos spin 
chain, following to a large extent the procedure of [7]. However, as these authors pointed 
out, we should stress that the underlying /i-expansion constitutes a somewhat ad- hoc pro- 
cedure and does not seem to shed much light on the fundamental origin of these conser- 
vation laws. One wonders whether the complicated limiting procedure is really necessary 
and whether these invariants could not be generated in a simpler way, from an intrinsic 
spin-chain formulation. In addition, we could ask whether explicit expressions for these 
Hamiltonian conservation laws could be written, in analogy with those of the XXX model 
[23] . We definitely see a similar pattern emerging but the expressions for the relative coef- 
ficients of the various terms appear rather complicated. Finally, a brute force computation 
of Hi and H3 seems to hint at the existence of an odd set of involutive charges which does 
not obey the yangian symmetry, and for which we still lack a generating function. 



Appendix A. A direct computation of [H 2 , H4] 

In this appendix, we show that the commutator [ H 2 , H4 ] vanishes by calculating it 
explicitly. For compactness, let us start by expressing the conservation laws in the form 

N 

H 2 = hijPij (A.l) 

N N ^ 

Ha = — ^2 hijkiPijki + ^2 fijPij ' (A-2) 

i,j,k,£=l *,j=l 



where 



fHj = -±- (A.3) 

hijki = (A.4) 

/« = Jr-|^-|(^-D^--|- (A*) 

ij ij ij 
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A direct calculation yields the following commutator 

JV 



[H2 , H4] =8 ^ (hij — hi m )hjke m Pijkem 

i,j,k,£,m 
N 

— 4 ^2 (hik - hje)hijkePijPk£ (A.6) 



i,j,k 
N f N 



JV r JV 

4^ < 2 ^ (/i^ - hu)hiju - {hik - hjk)fij \Pijk ■ 

i.i.k ^ 1=1 ' 



i,j,k y I 

£^i,j,k 



Defining now the cyclic sum operator 



cyclic 

^ ] f{ x il ■ ■ ■ x ik ) = f( x ii ■ ■ ■ x ik) "I" f{ x i2 > ^3 • • • ^ife ' ^1) 
{il...i fc } 

+ • • • + f( X ik > X ii • • • X ifc-l ) 



(A.7) 



and using the fact that the exchange operators in (A.6) are invariant under cyclic permu- 
tations of their indices, we can write 

N cyclic 



JV ^ cyclic ,. ^ 

[H 2 , H4] =8 ^ - ^ < (/ijj - h irn )hjk£m \PijUm 
i,j,k,e,m {i,j,k,£,m} ^ 

V ^ cyclic cyclic ^ ^ 
~ 4 4 Yl E| ^ ~ h 3i) h ijk£ >PijPk£ (A.8) 

n n U I) I A fl, <1 V J 



{M} 

jV 1 cyclic .. JV 



JV ^ ^y^HV> ^ IV N 

4^ - ^ <2 ^ (ha - hk£)hijk£ - (hik - hjk)fij?Pijk ■ 

i,3,k {i,j,k} ^ fcl ' 



This commutator will therefore vanish if one can prove that 

cyclic 



F 1 = ^ \ ( h ij ~ h im)hj k £m > = (A.9) 
{i,j,k,£,m} ^ * 
cyclic cyclic , 

F 2 =J2Y1 ~ h je) h ijk£j = (A.10) 



{k,£} 
cyclic r TV 



^3 = ] 2 ( hii ~ h k£)hijk£ + (hik - h jk )fij > = . (A.ll) 



£^i,j,k 
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The first condition (A. 9) is shown to be satisfied in the following manner. First, we extract 
a cyclic invariant from the sum 

cyclic , , _ N 

1 



cyclic 



y l ^_ , 



XijX jkXktXtmXmi r . . , , , [XijXmj %mi%mj 
{i,j,k,t,m} v 

cyclic ^ ^ 

, \ r I ^mi ^ij I 

= hijklm 2_^i 1 | 



{i,j,k,£,m} 

The next few steps are just basic algebra 

cyclic , o rJl 



f X ■ — X- ■ 1 

Fi = hijktm Yl \ ™ . r tJ . \ 

{i,j,k,e,m} ^ X *3 X m3X mt J 

Cyclic ( ( _i_ \ _ o ^ 

= /lyfc/m 2^ i ~ ~ ~ r ( A - 13 ^ 

{i,j,k,i,m} 1 x *j x mj x ™ ) 

cyclic , _ . 

\ I X m -f- Xj ZXi I 

- fiijkem 2^ i f • 



{i,j',fe^,m} 

Extracting once more the cyclic invariant hijktm gives 

cyclic , . 

^1 = hijktm Y 1 ( Xm + X 3 ~ 2 %i)XjkXktXtm > • (A. 14) 

{i,j,k,t,m} 

This cyclic sum can then be shown to vanish by plainly writing down all of its terms. 
The second condition (A. 10) can be proved to hold in a similar fashion. Establishing the 
vanishing of F 3 is a bit more tricky however. The main steps are as follows. First, we 
write F3 explicitly: 

cyclic f « r / 1 j ^ J J 



1 3 3 I 3 3 

x ij x jk I e=1 \ x kt x ti x H x tk/ X kj x ji X ij X jk- 

e&.h (A. 15) 

V4 ^/U 3 4 3 4 3 

We start by using the summation identity (B.17) in order to simplify the first two terms 
and notice that the last two terms in the last parenthesis do not contribute. The reduced 
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expression is: 



cyclic 

f»= E 

{i,3,k} 



1 X^, 1 Xj -|- Xfc 



X 



ki 



+ 



2 x^j 
2 



+ 



O I o 



2 37 4 37 4 



(A.16) 



™4 2 
X ij X jk 



+ 7T- 



2 37 ^ 37 j 



™4 ™2 q „2 —I Q ^,2 ™2 



The cyclic sum will also cancel the last two terms in square brackets with the subsequent 
two terms so that we are left with 



cyclic 

*■» = E 

{i,3,k} 
1 



1 X 2 H 37 /j; 2 37 ^ 37 jf 2 37 2 37 j 

3 XijXjkXhi 3 x^jX^fe 



3 XijX jk%k 



cyclic ^ 2 cyclic 

^ (Xi + Xfc) + --y— 2~ ^ XiXj(x jk — X ki ) 



{i,j,k} 



o i2 ™ 2 ™ 2 

5 ij jk ki {iJM 



2 (xi + + x fc ) 2 (x;x^ - xf Xj + Xjx\ - x^x k + x k x\ - xix\) 



+ 



9 9 9 
^ij^jk^ki 



3 Xjj XjfcX/jj 3 

2 (Xj -(- Xj -f- Xfc) 2 (Xj -f- Xj -\- X/ c )XjjX jk^k 



3 XjjXj^x^j 



+ 



/y»2 ™ 2 -yi 2 

x ij x jk x ki 



(A.17) 



We therefore see that #2 and H4 do commute, a fact which corroborates the validity of 
the dynamical ^-expansion used throughout this work. 



Appendix B. The zeroes of the Hermite polynomials: summation identities 

In this appendix, we briefly show how the lattice sites of the FP model, defined by 
(4.6), can be identified with the zeroes of the Hermite polynomial Hn(x) and then present 
a series of summation identities which are vital for the reduction of certain expressions. 
Following [14], consider then the Hermite differential equation 

H'^(x)-2xH' N (x) + 2NH N (x) = 0. (B.l) 

Letting Xj(z = 1 . . . N) denote the zeroes of H N (x) and evaluating (B.l) at an arbitrary 
zero xi gives 

H^(x £ ) = 2x e H' N (x e ) . (B.2) 
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Factorizing H N (x) in terms of its zeroes and substituting in (B.2) generates the identity 



N 1 

— = x j ■ 



(B.3) 



fc=l 



As already mentioned, a number of simple summation identities generalizing the previous 
one have already been discovered some time ago [22] . One can easily generate more com- 
plicated formulae. The general procedure is the following: to increment a power to the 
numerator of (B.3), we can proceed as follows 



N 



x, 



k=l 



%jk 



N 

E 



1 



N 



k=i X i k k=i X i k 



fc=l 

k^j 



N 



=► ]T^ = z 2 -(iV-l). 

fc=i Xifc 
Mi 

On the other hand, to increase a power in the denominator, the procedure is 



(B.4) 



N 



N 



N 



k=i x i k k=i Xik 



k^j 
N 



k^i 



k^i,j 

that is 



E— = -- 

fc—l X jk%ik %ij 



k^i,j 
N 



Xij 



El 1 _ y> Xifc — ^jfc 1 1 

^ _^ •^jk'L'ik ^ji >^ij fc— 1 jk^ki >t-ij 

k^i,j 



N N 

E — -£- 

* y . , # . 



fc=l Xjfc fc=l 



H 
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N 



, - 1 - 2 ' 

fc=1 X jk X ik %ij 



E 

fc=i 

k^i,j 



(B.5) 



(B.6) 



Using such procedures, one can generate a whole set of summation identities, the most 
useful of them being given by 



N 

E — 

, 1 %ik 



Xi 



k=i 

N 



Et = !("-')-5 



/.•— 1 3 



X 2 

3 1 



(B.7) 



(B.8) 
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k=i ik 
N i 



k=i Xik 
k^i 

N 



k=l 
k^i 

N 



Xik 



E 



k=i Xik 
k^i 

N o 

xi 



E^k _ 

k=l ^ 
k^i 

N 



1 3 



§("-l)-l] 



E — = 4-(Af-2)x, 



fc=l 
N 1 

E — 

N 

E -Ar 

fc=i XlfeX fci 

N 1 

E 



1 - 



1 2iV + 1 - x\ 



x 



X 



■I'J 



1 27V + 1 - xj 1 Xj 



k=i 



XlkX k j 



X 



4 

T4~ 



2 x^- 



(B.10) 



(B.11) 



(B.12) 



(B.13) 



(B.14) 



(B.15) 



(B.16) 



(B.17) 



We now notice that (B.9) is identical to (4.6), the condition for the FP model to be 
y[s«(n)]-symmetric. This proves that the lattice sites of the FP model with N spins are 
actually the zeroes of the Hn(x) polynomial. 
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